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(57) Abstract 

Apparatus and methods for detecting and correcting various fault conditions in an operating Coriolis effect mass flowmeter. The 
present invention receives information from an operating Coriolis mass flowmeter and compares the information to threshold signatures 
representing various fault conditions. When a fault condition is dectected, output signals are applied to inform an operator and to control 
the mass flow rate through the flowmeter to correct the fault condition. Specific methods of the present invention detect the presence of 
a crack in the flow tubes and stop the flow of material to prevent release of the material dirough cracked flow mbes. Other methods of 
the present invention detect the void fraction of material flowing through the flow tubes, compute a corrected actual mass flow rate, and 
control the mass flow rate through the flowmeter to compensate for the effects of the void fraction. 
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METHOD AND APPARATUS FOR FAULT DETECTION AND 
CORRECTION IN CORIOUS EFFECT MASS FLOWMETERS 

5 RELD OF THE INVENTION 

This invention relates to Goriolis effect mass flowmeters used to measure 
flow of a material through a pipe or conduit. In particular, this invention relates to 
detection and correction of various fault conditions which may arise in the operation 
of a Goriolis effect mass flowmeter. 

10 

PROBLEM 

It is known to use Goriolis effect mass flowmeters to measure mass flow and 
other information for materials flowing through a conduit. As disclosed in U.S. Pat. 
Nos. 4,491,025 issued to J. E. Smith, et a!., of January 1, 1985, and Re. 31,450 to 

15 J. E. Smith of February 11, 1982, these flowmeters have one or more flow tubes 
of a curved configuration. Goriolis effect mass flowmeters having straight flow 
tubes are also well known as shown in U.S. Pat. No. 4,109,524 Issued August 29, 
1978 also to J. E. Smith, et al. Each flow tube configuration in a Goriolis mass 
flowmeter has a set of natural vibration modes, which may be of a simple bending, 

20 torsional or coupled type. Each flow conduit is driven to oscillate at resonance in 
one of these natural modes. The natural vibration modes of the vibrating, material 
filled system are defined in part by tlie combined mass of the flow tubes and the 
material within the flow tubes. Material flows into the flowmeter from a connected 
conduit on the inlet side of the flowmeter, is directed through the flow tube or tubes, 

25 and exits the flowmeter to a conduit connected to its outlet side. 

When there is no flow through the flowmeter, all points along the flow tube 
oscillate with identical phase due to an applied driver force. As material begins to 
flow, Goriolis accelerations cause each point along the flow tut>e to have a different 
phase. The phase on the inlet side of the flow tube lags the driver, while the phase 

30 on the outlet side leads the driver. Sensors are placed on the flow tube to produce 
sinusoidal signals representative of the motion of the flow tube. The phase 
difference between two sensor signals is proportional to the mass flow rate of 
material through the flow tube. 
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A complicating factor In this measurement Is that the density of typical 
process material varies. Changes in density cause the frequencies of the natural 
vibration modes to vary. Since the flowmeter's control system maintains resonance, 
the oscillation frequency varies in response. Mass flow rate in this situation is 
5 proportional to the ratio of phase difference and oscillation frequency. 

The Coriolis flowmeter of U.S. Pat. No. Re. 31,450 avoids the need of 
measuring both phase difference and oscillation frequency. Phase difference is 
determined by measuring the time delay between level crossings of two sinusoidal 
output signals from the flowmeter. When this method is used, the variations in the 
10 oscillation frequency cancel, and mass flow rate is proportional to the measured 
time delay. This measurement method is frequently referred to as a time delay 
measurement. 

In many applications of Coriolis mass fbwmeters the material flowing through 
the flowmeter and flow tubes may be corrosive or under high pressure. A crack In 

1 5 the flow tubes of such a mass flowmeter could release corrosive material. In these 
applications it is important to minimize the possibility of leakage of the material. 
One method to minimize such risks is to encase the flowmeter so as to contain any 
material which may escape in case of a loss of pressure integrity in the flow tubes. 
This solution presents many additional problems depending on the materials used 

20 and the operating pressure levels. The complexity and cost of encasing flow tubes 
containing corrosive materials, or containing materials at very high pressures can 
be prohibitive in many applications of Coriolis effect mass flowmeters. 

Rather than attempt to contain the flow tubes in case of such a fault, it is 
preferred to detect the fault condition before the fluid can be released to the 

25 surrounding environment. Eariy detection of such a fault permits the flow of the 
material to be stopped before leakage occurs. Others have tried to address this 
problem by monitoring the sound waves emanating from the vibrating tubes 
(acoustic emission spectroscopy) or by measurement of changes in resistivity of the 
flow tubes (electronic fingerprint method). Both of these methods can detect cracks 

30 in the flow tubes, however they are similariy sensitive to flow tube surface 
properties such as conx)sion. To distinguish the growth of a crack from other 
surface conditions by use of these methods requires sensitive electronic 
measurement and skilled interpretation of the measurement results. The complexity 
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of these measurements has rendered these methods of little use In practical 
commercial applications of flowmeters. 

Another problem with prior designs arises when air, or other foreign material, 
infiltrates the material flowing through a Coriolis effect mass flowmeter. These 
5 materials alter the vibrational characteristics of the flowmeter flow tubes in such a 
manner that the computed mass flow is erroneous. Prior designs neither detect nor 
correct these errors in mass flow measurements. 



S0|,MTION 

1 0 The above and other problems are solved by the present invention which 

comprises methods and apparatus for detecting various fault conditions in an 
operating Coriolis effect mass flowmeter as well as methods and apparatus for 
correcting those detected faults. The present invention provides for the detection 
of cracks in the flow tubes early enough to permit the flow of the material to l>e 

15 stopped before a release of the flowing material to the surrounding environment. 
In addition, the present invention permits early detection and correction of other 
conditions relating to the flowing material which may indicate other types of system 
failure. Changes in computed mass flow, material density, void fraction, or mass 
fraction beyond threshold limits may represent system failures in certain applications 

20 of mass flowmeters. The present invention comprises methods and apparatus for 
monitoring the operation of Coriolis effect mass flowmeters to permit early detection 
of a crack in the flow tubes or the aforementioned system failures. Additionally, the 
present invention comprises methods and apparatus for controlling the material flow 
within the flowmeter to correct the detected faults or to stop the flow of material to 

25 prevent harmful release of the flowing material. 

In a normally operating flowmeter, under static conditions, the drive power 
required to vibrate the flow tubes and the resulting, resonant vibration frequency of 
the flow tubes are constant. Changes in the required drive power or resonant 
frequency of the vibrating flow tubes indicate a change in the conditions in which 

30 the fbwmeter is operating. The changes in the drive power and resonant frequency 
are analyzed by the methods and apparatus of the present invention to detect a 
crack in the flow tubes or other operational anomalies. Each fault detectable by the 
present invention exhibits a characteristic "signature" in the changes over time of 
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drive power and/or resonant frequency. For example, a change in the mass flow 
rate of material in the flow tubes can cause a change of the resonant frequency of 
the flow tubes without a significant change in the required drive power. However, 
a sudden decrease in resonant frequency of the flow tubes when coupled with an 
5 Increase in required drive power is indicative of a crack in the flow tubes. Other 
changes in the operating conditions exhibit similariy recognizable "signatures" in the 
operating parameters of the flowmeter and are detectable by the present invention 
as operating faults or anomalies. The "signatures" are defined in terms of threshold 
values for the drive power and frequency as well as the rate and cun/atui« of 
1 0 changes over time (first and second derivatives with respect to time). 

The present invention comprises apparatus for sensing certain operational 
parameters of a Coriolis effect mass flowmeter including; the drive power required 
to vibrate the flow tubes and the resonant frequency at which the tubes vibrate in 
response to the applied drive power. The apparatus of the present invention further 
1 5 comprises output signals for intorming an operator of the occun-ence of a particular 
fault condition in the operation of the flowmeter as well as output signals for 
automatically adjusting or stopping the flow of material through the flowmeter. 
These latter output signals are applied to control pumps and valves associated with 
the flow of material in the flowmeter. 
20 The present invention further comprises methods for detennlning changes 

in the sensed operational parameters over time. These methods determine the 
instantaneous values of certain sensed operational parameters as well as the rates 
of change, (slope and curvature), of each sensed operational parameter. These 
methods analyze the present values and changes of the sensed operational 
25 parameters to determine that a particular fault condition is occurring. Each fault 
condition detectable by the present invention is recognized by a unique "signature" 
in the values and changes of the sensed operational parameters. Signature values 
for the operational parameters which indicate various fault conditions are stored 
within the apparatus of the present invention. To the extent that the values and 
30 changes of the monitored operational parameters match these stored signatures, 
the aforementioned output signals are generated by the apparatus of the present 
invention to indicate the occurrence of the related fault condition or to control the 
material flow through the flowmeter. 
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in particular, the present invention detects the presence of a crack in the flow 
tubes of a Coriolis flowmeter. The presence of a cradk is recognized by the 
"signature" decrease in resonant frequency of the vibrating flow tubes coupled with 
increasing drive power required to vibrate the flow tubes. This early detection of 
5 cracl( faults permits an operator or control circuit to stop the flow or reduce the 
pressure in the flow tubes so as to prevent hannful leakage of the material flowing 
in the flowmeter. 

The present invention also recognizes operating conditions which may 
represent faults in some flowmeter applications such as mass flow changes, 

1 0 material density changes, void fraction changes, and mass fraction changes. For 
example, an increase in resonant frequency coupled with an increase in drive power 
indicates an increase in void fraction of the material flowing in the flowmeter. In 
some applications, this may be considered a fault condition indicating that, for 
example, excess air bubbles or pockets are fonnlng in the flowing material or that 

1 5 the output of the sensor is in enor. Output signals provided by the apparatus of the 
present invention are generated to wam an operator, to automatically stop the ftow, 
or to automatically adjust the mixture of the material flowing in the conduit. 
Additional output signals are applied to automatically adjust the flow rate of the 
material to minimize the error for increased void fraction. 

20 In addition to detection of fault conditions, the methods and apparatus of the 

present invention automatically correct the measured flow rate In the operating 
flowmeter to account for void fraction or mass fraction in the flowing material. 
Information stored within the apparatus of the present invention is used to compute 
the void fraction as a function of the measured values and computed rates of 

25 change in vibrating frequency and drive power. TTiis comFXJted void fraction 
coupled with the measured ftow rate is then used to compute an enx>r percentage. 
The error percentage is eipplied to the measured mass flow rate to compute a 
corrected mass flow rate which is used for operator display or automatic flow 
control as discussed above. 

30 In sum, the present invention provides methods and apparatus for eariy 

detection of various fault conditions in the operation of a Coriolis effect mass 
flowmeter. The present invention provides fault detection and correction with 
improved reliability and with less complexity than previous methods or apparatus. 
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The present invention also provides methods and apparatus for automatically 
correcting the mass flow rate based on measured flow rate and computed void 
fraction. 

5 BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 depicts a typical Corioiis effect mass flowmeter with the fault detection 
apparatus of the present invention; 

FIG. 1 A depicts the Corioiis effect mass flowmeter of FIG. 1 connected to a 
controllable pump and controllable valve in a typical closed loop control application 
10 configuration; 

FIG. 2 shows additional detail of the fault detection instrumentation of the 
present invention; 

FIG. 3 is a graph which depicts the characteristic signature of the frequency 
and drive power curves as a function of time with a crack in the flow tubes; 
15 FIG. 4 is a graph which depicts the characteristic signature of the frequency 

and drive power curves as a function of mass flow rate within the flow tubes; 

FIG. 5 is a graph which depicts the characteristic signature of the frequency 
and drive power curves as a function of material density flowing within the flow 
tubes; 

20 FIG. 6 is a graph which depicts the characteristic signature of the frequency 

and drive power curves as a function of void fraction or slug flow of the material 
flowing within the flow tubes; 

FIG. 7 is a graph which depicts the characteristic signature of the frequency 
and drive power curves as a function of mass fraction of the material flowing within 
25 the flow tubes; 

FIG. 8 is a graph which depicts the eror percentage to be applied to 
measured mass flow rate to compute corrected mass flow rate; 

FIG. 9 is a flowchart depicting an overview of the steps comprising the 
method of the present invention; 
30 FIG. 10 is a flowchart depicting additional detail of the step of FIG. 9 for 

sampling operational parameters; 

FIG. 11 is a flowchart depicting additional detail of the step of FIG. 9 for 
storing sampled operational parameters; 
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FiG. 12 is a flowchart depicting additional detail of the step of FIG. 9 for 
computing slopes and curvatures of changes in stored operational parameters; 

FIG. 13 Is a flowchart depicting additional detail of the step of FIG. 9 which 
tests for a cracked flow tubes and stops material flow accordingly; 
5 FIG. 14 is a flowchart depicting additional detail of the step of FIG. 9 which 

computes void fraction corrected mass flow rate and adjusts material flow 
accordingly; and 

FIG. 15 is a flowchart depicting additional detail of the step of FIG. 9 which 
tests for other operating fault conditions and controls mass flow rate accordingly. 

10 

DETAILED DESCRIPTION OF THE DRAWINGS 

Apparatus: 

A typical Coriolis effect mass flowmeter 10 is illustrated in FIG. 1 as having 
two cantilever mounted flow tubes 12, 14 affixed to a manifold body 30 so as to 

1 5 have substantially identical spring constants and moments of inertia about their 
respective out-of-phase bending axes W-W and W'-W. 

A drive coil and magnet 20 are mounted at a midpoint region between the 
top portion 130 and 130' of flow tubes 12, 14 to oscillate flow tubes 12, 14 out of 
phase about axes W-W and W'-W'. Left sensor 16 and right sensor 18 are 

20 mounted near the respective ends of the top portions of flow tubes 12, 1 4 to sense 
the relative movement of flow tubes 12, 14. This sensing may be done in many 
ways including by measuring the movement of the top ends of the flow tubes 12, 
14 through their zero crossings or some other pre-defined point. Flow tubes 12 
and 14 have left side legs 131 and 131 ' and right side legs 134 and 134'. The 

25 side legs converge downwardly toward each other and are affixed to surfaces 120 
and 120' of manifold elements 121 and 121'. Brace bars 140R and 140L are 
brazed to the legs of flow tubes 12, 14 and serve to define the axes W-W and W'- 
W' about which the flow tubes oscillate out of phase when driver 20 is energized 
over path 156. The position of axes W-W and W-W' is determined by the 

30 placement of brace bars 140R and 140L on flow tube side legs 131 , 131 ' and 134, 
134'. 

Temperature detector 22 is mounted on side leg 131 of flow tube 14 to 
measure the flow tube's temperature and the approximate temperature of the 
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material flovwng therein. This temperature information is used to determine 
changes in the spring constant of the flow tubes. Driver 20. sensors 16 and 18 and 
temperature detector 22 are connected to mass flow Instrumentation 24 by 
paths 157, 158 and 159. Mass flow instrumentation 24 may include a micro- 
5 processor which processes the signals received from sensors 16, 18, and 22 to 
determine the mass flow rate of the material flowing through flowmeter 10 as well 
as other measurements, such as material density and temperature. Mass flow 
instaimentation 24 also applies a drive signal over path 156 to driver 20 to oscillate 
tubes 12 and 14 out-of-phase about axes W-W and W'-W. 

10 Manifold body 30 is formed of casting 150, 150'. Casting elements 150, 

150' are attachable to a supply conduit and exit conduit (not shown), by flanges 
103, 103'. Manifold body 30 diverts the material flow from the supply conduit into 
flow tubes 12, 14 and then back into an exit conduit. When manifold flanges 103 
and 103' are connected via inlet end 104 and outlet end 104' to a conduit system 

15 (not shown), carrying the process material to be measured, the material enters 
manifold body 30 and manifold element 110 through inlet orifice 101 in flange 103 
and is connected by a channel (not shown) having a gradually changing 
cross-section in casting element 150 to flow tubes 12, 14. The material is divided 
and routed by manifold element 121 to the left legs 131 and 131 ' of flow tubes 14 

20 and 12, respectively. The material then flows through the top tubes elements 130, 
130' and through the right side legs 134 and 134' and is recombined into a single 
stream within flow tube manifold element 121 '. The fluid is thereafter routed to a 
channel (not shown) in exit casting element 150' and then to exit manifold element 
110'. Exit end 104' is connected by flange 103' having bolt holes 102' to the 

25 conduit system (not shown). The material exits through outlet orifice 101* to return 
to the flow in the conduit system (not shown). 

Mass flow instrumentation 24 analyzes signals received on paths 157, 158, 
and 159 and generates standard output signals on path 155 to indicate mass flow 
rates utilized by a control system or operator for monitoring and control of the mass 

30 flow rate through the associated conduit system (not shown). In addition, mass ftow 
instrumentation 24 generates output signals on paths 160, 161, 162. and 163 
indicating present resonant frequency of the flow tubes, present drive power 
required to drive the flow tubes, present temperature of the flow tubes, and 
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computed mass flow rate, respectively. These signals are received by fault 
detection Instmmentation 26 and analyzed to determine the occurrence of faults in 
the operation of flowmeter 10. These signals (and others detectable in the 
operation of flowmeter 10) are hereinafter referred to collectively and individually as 
5 "operating parameters". 

Mass flow instrumentation 24 comprises standard apparatus and methods 
known to the art. A Micromotion Remote Flow Transmitter model number 9739 is 
exemplary of such a device. Although a signal indicative drive power is not 
normally provided as an output signal from the model 9739, one skilled in the art 

10 would readily recognize that such a signal is available from the 9739 instrument. 
A frequency signal is commonly available from the model 9739 since the design is 
such that fluid density is measured. An electronics technician could easily modify 
such a device to apply appropriate signals indteative of flow tube vibration frequency 
and required drive power to conductore directed out of the instrument connecting 

15 to a fault detection instrumentation 26. 

Upon detection of fault conditions, fault detection instrumentation 26 applies 
signals to pump 180 over path 170 and to valve 181 over path 171 to correct the 
mass flow rate through flowmeter 1 0 or to stop the material flow to prevent release 
of the material. Pump 180 is controlled by fault detection instrumentation 26 to 

20 increase or decrease the pressure at inlet orifice 1 01 of flowmeter 1 0 to connect a 
detected fault. Valve 181 is controlled by fault detection instrumentation 26 to 
increase or decrease backpressure at outlet orifice lOV of flowmeter 10 to correct 
a detected fault. Fault detection instrumentation 26 also applies signals on path 
172 to inform the operator of the occurrence of a fault condition detected in the 

25 operation of flowmeter 10. The signals applied to path 172 are messages displayed 
on panel 182 to inform the operator of the need for manual corrective action or to 
inform the operator as to the status of corrective control performed by the fault 
detection instrumentation 26. 

For example, fault detection instrumentation 26 determines whether a crack 

30 is present in the flow tubes by analyzing the resonant frequency and drive power 
infomiation received on paths 160 and 161 respectively. If a crack is present in the 
flow tubes, fault detection Instrumentation 26 generates an output signal on path 
170 to decrease pressure generated by pump 180 at inlet orifice 101 of flowmeter 
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1 0. This slows or stops the flow of material through the flow tubes before the crack 
results in a loss of pressure integrity thus preventing release of the material from 
within the flow tubes 12 and 14. 

FIG. 1A depicts flowmeter 10 and fault detection instrumentation 26 with 
5 pump 180 and valve 181 connected as they would be in a typical application of 
flowmeter 10. Pump 180 draws material 192 from source 191 and delivers it to the 
Inlet orifice 101 of flowmeter 10 through conduit 190. The inlet pressure created 
by pump 180 is controlled by fault detection instrumentation 26 by applying signals 
over path 170 to increase or decrease the speed of pump 180. Material delivered 
10 to inlet orifice 101 of flowmeter 10 moves through the meter to outlet orifice 101 ' 
into conduit 190 and finally through valve 181 to its destination. The backpressure 
created by valve 181 is controlled by fault detection instrumentation 26 by applying 
signals over path 171 to open or close valve 181. 

It will be recognized by those skilled in the art that the methods and 
15 apparatus of the present invention may be advantageously applied to any 
configuration of flowmeter. The particular flowmeter configuration discussed above 
with respect to FIGS 1 and 1A is exemplary of all such flowmeter configurations. 
In particular, straight flow tube, U-shaped flow tube, in-egular and other 
configurations of flowmeters may be combined with the methods and apparatus of 
20 the present invention for detecting and con-ecting operating fault conditions. 

FIG. 2 shows additional detail of the components within fault detection 
instoimentation 26. Analog signals applied to paths 160 through 163 by mass flow 
instrumentation 24 of FIG. 1 are converted to digital values by analog to digital 
(A/D) converters 203 through 206. Signals representing the converted digital values 
25 are applied to paths 220 through 223 for processing by CPU 200. The signals 
applied to paths 220, 221, 222, and 223 are representative of present resonant 
frequency of the vibrating flow tubes, present drive power required to vibrate the 
flow tubes, present temperature of the flow tubes, and computed mass flow, 
respectively. 

30 CPU 200 periodically samples the digital signals representing present 

frequency, drive power, temperature, and computed flow rate from paths 220 
through 223 respectively. At each periodic sampling, CPU 200 averages multiple 
readings of the digital values applied to paths 220 through 223 and writes the 
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averaged values in buffers In RAM memory 201 over path 229. This averaging aids 
in filtering analog signal noise which may be present on paths 160 through 163. 
Averaged values for frequency, drive power, temperature, and computed flow rate 
are stored In buffers In RAM memory 201 in a circular manner so that when the 
5 buffers are filled the oldest previously stored values are overwritten by the newest 
values to be stored. The circular buffers contain the values of the stored operating 
parameters over a period of time. 

CPU 200 computes the first and second derivatives with respect to time for 
each stored operating parameter to determine the slope and curvature of any 
0 changes over time in the operation of flowmeter 10 of FIG. 1. CPU 200 reads 
stored averaged values from RAM memory 201 over path 229. The difference 
between two averages stored sequentially In time is the slope of the changes In that 
operating parameter (first derivative of the sampled operating parameter). An 
example of this conputation is described wherein fo through f^ are presumed to be 
5 the stored averages for frequency of flow tube vibrations over N+1 periods of 
reading by CPU 200 and wherein S, through are the computed slopes for the 
changes in the values of fo through f^ for each sampling period. 
S, = (f, - fo) 

Or more generally: S, = (f, - f,.,) where x - 1 through N. 

The difference between two such computed slopes indicates the curvature 
of the operating parameter (second derivative of the sampled operating parameter). 
An example of this computation is as follows wherein S, through are presumed 
to be the slopes of the changes in periodically stored averages as discussed above 
and wherein through C^ are the computed curvatures for the changes In the 
slopes S, through S^. 
Ca = (Sa - S,) 

Or more generally: C, = (S, - S,.,) where x - 2 through N. 

ROM memory 202 contains the threshold "signatures" for each fault 
detectable by the present invention. The Information stored In the "signatures" 
relates to threshold values for present values of operational parameters as well as 
threshold values for slope and curvature of changes in operational parameters. The 
infonnation is recorded for each operational fault detectable by the apparatus and 
methods of the present invention. The "signature" thresholds are unique to the 
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particular characteristics of the flowmeter, the nominal physical properties of the 
material flowing through the flowmeter, and the nominal flow rates of the material 
through the flowmeter. These threshold values are predetermined by the user or 
manufacturer of the flowmeter and programmed into ROM memory 202. In an 
5 altemative exemplary embodiment, CPU 200 may copy these values from ROM 
memory 202 into RAM memory 201 and interact with the operator to modify the 
threshold "signatures" for a particular customized need in the application of the 
flowmeter. 

ROM memory 202 contains additional information predetermined at time of 
10 manufacture or installation of flowmeter 10 pertaining to correction of measured 
mass flow rate to account for detected changes in material void fraction. As 
discussed below, increases in void fraction create errors in the measurement of 
mass flow rate. The error factor (error percentage) may be detemiined by methods 
of the present invention utilizing "signature" information which encodes the error 
1 5 percentage as a function of measured mass flow rate and computed void fraction. 
The "signature" information representing the relationships between measured mass 
flow rate, computed void fraction, and resultant en^or percentage Is unique to the 
particular flowmeter 10 and the particular application. This information is stored in 
ROM memory 202 with the fault detection signatures. The methods of the present 
20 invention which utilize this information to correct measured mass flow rate are 
discussed below. 

CPU 200 compares the stored average values of operating parameters, the 
computed slopes, and the computed curvatures to corresponding threshold values 
read over path 228 from the "signatures" stored in ROM memory 202. The 

25 signature information indicates specific threshold values for the operating 
parameters, the computed slopes, and the computed curvatures which indicate the 
occun^ence of particular fault conditions in the operation of flowmeter 10 of FIG. 1 . 
CPU 200 compares the signature threshold Information read from ROM memory 
202 to the con^esponding averaged values in RAM for frequency, drive power, 

30 temperature, and measured mass flow rate. CPU 200 also compares the signature 
threshold values in ROM memory 202 to the computed slopes and curvatures for 
each parameter. These comparisons determines whether each of various fault 
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conclitions has ocxnjrrecl In the operation of flowmeter 10 of FIG. 1. The methods 
applied to make this comparison are discussed below in further detail. 

CPU 200 applies signals to path 170 to control the speed of pump 180 and 
applies signals to path 171 to open or close valve 181 In response to detection of 
5 a fault condition. In addition, CPU 200 displays messages on panel 182 over path 
172 to inform the operator regarding controls and status of the operation of 
flowmeter 10. As discussed below, when CPU 200 detects a crack in the flow 
tubes 12 and 14 of flowmeter 10, it controls pump 180 and valve 181 to stop the 
flow of material and relieve the pressure within flowmeter 10. Similarly, CPU 200 

1 0 may modify the flow rate to compensate for other fault conditions such as due to: 
faulty mass fk)w, faulty material density, faulty void fraction, or faulty mass fraction. 
Each of these fault conditions is sensed by CPU 200 by comparison of the 
signature threshold infomiation stored in ROM memory 202 with present parameter 
values stored in RAM and computed sbpes and curvatures for each parameter as 

1 5 discussed above. Each fault condition is defined by a set of threshokJ values stored 
in ROM memory 202. The threshold values relate to the present values of the 
several operating parameters as well as the slope and curvature of changes in the 
several operating parameters. 

20 Sionature Information Samples: 

FIG. 3 is a graph depicting typical changes in resonant frequency of the 
vibrating flow tubes and changes in drive power required to vibrate the flow tubes 
as a crack forms and grows through the wall of a flow tube. The abscissa of the 
graph in FIG. 3 is time in seconds. The curve 300 labeled "Pressure" has it's 

25 ordinate scale labelled on the left side of the graph. This cun/e 300 represents the 
pressure within the flow tubes over time as the crack in the flow tubes grows over 
time. As can be seen in the graph, pressure wrthtn the flow tubes remains constant 
until time 180 seconds shown by a dashed line lat>elled 312. Pressure begins to 
drop at this point due to release of the flow material through the crack in the flow 

30 tube walls. 

The curve 301 labelled "Drive Power^ also has it's ordinate scale labelled on 
the left side of the graph. Curve 301 represents the drive power required to vibrate 
the flow tubes at the resonant frequency. The curve 302 labelled "Frequency" has 
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it's ordinate scale labelled on the right side of the graph. Curve 302 represents the 
change In the resonant frequency of the vibrating flow tubes. The drive power and 
frequency remain relatively constant as the crack begins to grow through 
approximately time 80 seconds shown by a dashed line labelled 310. Frequency 
5 curve 302 begins to decrease and drive power curve 301 begins to increase at this 
time. At time 80 seconds labelled 31 0, CPU 200 may determine that a crack is 
growing in the flow tubes by matching the frequency decrease and drive power 
Increase with appropriate signature information stored in ROM memory 202. A non- 
linear increase in the required drive power as exhibited in curve 301 and a non- 
10 linear decrease in the resonant frequency as exhibited In curve 302 are signature 
characteristics of a growing crack In the flow tubes of a flowmeter such as 
flowmeter 10 of FIG. 1 . As can be seen by the relationship of the curves 300, 301 , 
and 302, of FIG. 3, methods of the present invention can detect the fonnation of a 
crack In flow tubes 12 or 14 as soon as time 80 seconds (labelled 310), well before 
15 the release of pressure from within flow tubes 12 or 14 as indicated by the line 
labelled 312 at time 180 seconds. 

The signature change In operating parameters of a forming crack fault 
depicted in FIG. 3 and discussed above Is clearly distinguishable from a normally 
operating flowmeter responding to changes In the mass flow of material (as 
20 depicted in FIG. 4). FIG. 4 depicts the signature change in operating parameters 
of a flowmeter normally responding to changes in material mass flow through the 
flowmeter. The abscissa of the graph in FIG. 4 is mass flow rate of the material 
flowing through the flow tubes of a flowmeter such as flowmeter 10 of FIG. 1 . The 
curve 401 labelled "Drive Power has it's ordinate scale labelled on the left side of 
25 the graph. This curve 401 represents the drive power required to vibrate the flow 
tubes at the resonant frequency. The curve 402 lat>elled "Frequency" has Ifs 
ordinate scale labelled on the right side of the graph. This curve 402 represents 
the change in the resonant frequency of the vibrating flow tubes. As shown by the 
curves of FIG. 4, in a normally operating flowmeter, drive power remains constant 
30 while mass flow rate increases but resonant frequency decreases while mass flow 
rate increases. Constant drive power as exhibited in curve 401 and a non-linear 
decrease in the resonant frequency as exhibited in curve 402 are signature 
characteristics of an increasing flow rate in the flow tubes of a flowmeter such as 
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fiowmeter 1 0 of FIG. 1 . Conversely, a constant drive power coupled with a non- 
linear Increase in resonant frequency indicates a decreasing flow rate in the flow 
tubes of a flowmeter 

FIG. 5 depicts the signature change in operating parameters of a fbwmeter 
5 normally responding to changes in density of the material flowing through the 
flowmeter. The abscissa of the graph in FIG. 5 is the density of the material flowing 
through the flow tubes of a flowmeter such as flowmeter 10 of FIG. 1. The curve 
501 labelled "Drive Power" has it's ordinate scale labelled on the left side of the 
graph. This curve 501 represents the drive power required to vibrate the flow tubes 

10 at the resonant frequency. The curve 502 labelled "Frequency" has irs ordinate 
scale labelled on the right side of the graph. This curve 502 represents the change 
in the resonant frequency of the vibrating flow tut>es. As shown by the curves of 
FIG. 5, drive power remains constant while material density increases but resonant 
frequency decreases while material density increases. Constant drive power as 

15 exhibited in curve 501 and a non-linear decrease in the resonant frequency as 
exhibited in curve 502 are signature characteristics of increasing material density 
in the flow tubes of a flowmeter such as flowmeter 10 of FIG. 1. Conversely, a 
constant drive power coupled with a non-linear increase in resonant frequency 
indicates a decreasing material density in the flow tubes of a flowmeter. 

20 The signature of changing material density as depicted in FIG. 5 is clearly 

distinguishable from a crack fault discussed above with reference to FIG. 3. The 
constant drive power exhibited by curve 501 in FIG. 5 Is clearly distinguishable by 
CPU 200 from the positively sloped, increasing drive power of curve 301 of FIG- 3. 
To distinguish the material density signature of FIG. 5 from the flow rate signature 

25 of FIG. 4, CPU 200 detects the negative curvature (second derivative) of frequency 
curve 402 in FIG. 4 from the positive curvature of frequency curve 502 in FIG. 5. 
The fault signature information stored in ROM memory 202 for a flow rate fault 
condition versus a material density fault condition indicates that CPU 200 must 
compare the curvature (second derivative) of the frequency operating parameter to 

30 determine which of the faults has occurred. 

FIG. 6 depicts the signature change in operating parameters of a flowmeter 
normally responding to changes in void fraction of the material flowing through the 
flowmeter. Void fraction is a measure of the air or gas content of the material 
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flowing through the flowmeter. The abscissa of the graph in FIG. 6 is the void 
fraction of the material flowing through the flow tubes of a flowmeter such as 
flowmeter 10 of FIG. 1. The curve 601 labelled "Drive Power has It's ordinate 
scale labelled on the left side of the graph. This curve 601 represents the drive 
5 power required to vibrate the flow tubes at their resonant frequency. The curve 602 
labelled "Frequency" has it's ordinate scale labelled on the right side of the graph. 
This curve 602 represents the change in the resonant frequency of the vibrating 
flow tubes. As shown by the curves of FIG. 6, both drive power and resonant 
frequency Increase while void fraction increases. Increasing drive power as 
10 exhibited in curve 601 and increasing resonant frequency as exhibited in curve 602 
are signature characteristics of increasing void fraction in the flow tubes of a 
flowmeter such as flowmeter 10 of FIG. 1. Conversely, decreases in both drive 
power and resonant frequency indicates a decreasing void fraction in the flow tubes 
of a flowmeter. 

1 5 FIG. 7 depicts the signature change in operating parameters of a fbwmeter 

normally responding to changes in mass fraction of the material flowing through the 
flowmeter. Mass fraction is a measure of the heavier material content in the 
material flowing through the flowmeter. The abscissa of the graph in FIG. 7 is the 
mass fraction of the material flowing through the flow tubes of a flowmeter such as 

20 flowmeter 10 of FIG. 1. The curve 701 labelled "Drive Power has it's ordinate 
scale labelled on the left side of the graph. This curve 701 represents the drive 
power required to vibrate the flow tubes at the resonant frequency. The curve 702 
labelled "Frequency" has it's ordinate scale labelled on the right side of the graph. 
This curve 702 represents the change in the resonant frequency of the vibrating 

25 flow tubes. As shown by the curves of FIG. 7, drive power increases while mass 
fraction increases and resonant frequency decreases while mass fractfon increases. 
Linearly increasing drive power as exhibited in curve 701 and linear decreasing 
resonant frequency as exhibited In curve 702 are signature characteristics of 
increasing mass fraction in the flow tubes of a flowmeter such as flowmeter 10 of 

30 FIG. 1 . Conversely, linear decreases in drive power and linear increases resonant 
frequency indicates a decreasing mass fraction in the flow tubes of a flowmeter. 

To distinguish the changing mass fraction signature of FIG. 7 from the crack 
fault signature of FIG. 3, CPU 200 must detect different curvature (second 
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derivative) of the frequency curves 302 in FIG. 3 and 702 In FIG. 7. Alternatively. 
CPU 200 must detect different curvature (second derivative) of the drive power 
curves 301 in FIG. 3 and 701 In FIG. 7. The fault signature information stored in 
ROM memory 202 for a crack fault condition versus a mass fraction fault condition 
5 indicates that CPU 200 must compare the curvature (second derivative) of the 
frequency or drive power operating parameter to determine which of the two faults 
has occurred. 

The precise "signature" threshold values, slopes, and curvatures for the 
operational parameters which indicate a particular fault are unique to each particular 

1 0 flowmeter and the nominal material composition, flow rate, and pressures within the 
flowmeter. The particular signature values are determined by testing at time of 
manufacture and stored in ROM memory 202. One skilled in the art will readily 
recognize that the signature threshold values may be copied to RAM memory 201 
by CPU 200 and then modified to customize for a particular application of the 

1 5 flowmeter. 

The ordinate scales of FIGS. 3 through 7 are expressed in percent change 
in the operational parameter. One skilled in the art will recognize that the threshokJ 
information may be equivalently expressed either as a percent change or as 
absolute threshold values. The methods of the present invention discussed below 

20 are expressed in terms of comparisons to absolute threshold values. 

FIG. 8 is a graph of information used in conjunction with information depicted 
in FIG. 6 to compute a corrected mass flow rate. The accuracy of the mass flow 
rate measured by mass flow Instrumentation is adversely affected by increases in 
void fraction of the material flowing through flowmeter 10. Void fraction, as 

25 discussed above, indicates the content of lighter material, such as air, mixed with 
the material flowing through flowmeter 10. FIG. 6 shows the signature infomnation 
which identifies an increasing void fraction as a function of frequency and drive 
power, spedfically: linear rising frequency coupled with a positively curved increase 
in drive power. 

30 In addition to other information discussed below, the signature information 

stored in ROM memory 202 includes an encoding of void fraction information 
depicted in the graph of FIG. 6. The void fraction information is used by the 
methods of the present invention to determine the current void fraction percentage 
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as a function of measured frequency and drive power. The void fraction tnfonnation 
of FIG. 6 is encoded into ROM memory, with the signature information, as 
coefficients of polynomial equations describing the curves of FIG. 6 or in a tabular 
form used to interpolate void fraction from frequency or drive power values. As 
5 described below with respect to FIG. 14, the void fraction information stored ROM 
memory 202 is used to compute the void fraction from frequency and drive power 
measurements. One skilled in the art will recognize that a tabular encoding of the 
functions represented In FIG. 6 Is preferred for use with low cost microprocessors 
which may possess limited mathematical computational power. 

10 Information shown in the graph of FIG. 8 is used to compute an error 

connection factor (error percentage) to be applied to the measured mass flow rate 
for determining a corrected mass flow rate. FIG. 8 shows six exemplary graphs of 
en-or percentage as a function of measured mass flow rate. Each of the six curves 
is labeled to indicate a con^esponding computed void fraction determined, as 

15 discussed above, utilizing information shown in FIG. 6 and stored with signature 
information in ROM memory 202. For example, if void fraction is computed as 
above to be 1.5%, the curve on FIG. 8 labeled 1.5 is used to determine the error 
correction factor (error percentage) to be applied given a mass flow rate. Using the 
curve labeled 1 .5 and the measured mass flow rate, the error percentage can be 

20 determined as shown in FIG. 8. 

in addition to other information discussed below, the signature information 
stored in ROM memory 202 includes an encoding of the error percentage 
information depicted in the graph of FIG. 8. The error percentage information is 
used by the methods of the present invention to determine the cun-ent error 

25 percentage as a function of measured flow rate and computed void fraction. The 
error percentage information of FIG. 8 is encoded into ROM memory 202, with the 
signature information, as coefficients of polynomial equations describing the curves 
of FIG. 8 or in a tabular form used to interpolate error percentage from measured 
mass flow and computed void fraction values. As described below with respect to 

30 FIG. 14, the en^or percentage infonmation stored in ROM memory 202 Is used to 
compute the corrected mass flow. One skilled in the art will recognize that a 
tabular encoding of the functions represented in FIG. 8 is prefen^ed for use with low 
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cost microprocessors which may possess limited mathematical computational 
power. 

Fault Signature In formation Format: 
5 Signature threshold Information is stored in ROM memory 202 of FIG. 2 for 

use by CPU 200 in comparing the cun-ent operational parameter values, slopes, 
and curves. The threshold information format is most easily understood when 
expressed as a data structure in a high level programming language such as "C". 
It is to be expressly understood that this expression is but one exemplary 

1 0 embodiment of the signature threshold information in ROM memory 202. Many 
equivalent forms of representing the reference threshold infonnation are within the 
spirit and scope of the present invention. The threshold information expresses the 
nominal values of a properly operating flowmeter by defining the upper an lower 
limits of the several operational parameters, of the slopes of the operational 

1 5 parameters, and of the curvatures of the slopes. The limits are encoded for each 
of several faults conditions detectable by the methods of the present invention. If 
the current values of the parameters, their slopes, and their curvatures, are within 
the defined limits then the corresponding fault condition has not arisen. 

struct reference_data { 

20 float NOMINALP; T nominal drive power for application V 

float NOMINALF; T nominal frequency for application V 

float MAXM; T maximum corrected mass flow V 

float MINM; /* minimum corrected mass flow V 
struct crack_deteclionJnfo { 

25 r 

* Cracks are indicated by positively curved increasing drive power 

* coupled with negatively curved decreasing frequency. 
*/ 

float MAXP; T maximum drive power V 

30 float MINF; /* minimum frequency */ 

float MAXPSLOPE; /* maximum drive power slope */ 

float MINFSLOPE; /* minimum frequency slope V 

float MAXPCURV; /* maximum drive power curvature V 

float MINFCURV; /* minimum frequency curvature V 
35 } CRACK; 

struct flow_correction Jnfo { 

r 

* Table describing signature curves of FIG. 6 

V 

40 float VOID_FRACTION[ ]; r table to compute void fraction V 
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r 

* Table describing plot of FIG. 8 
V 

float ERROR.PERCENTT 1; /* table to compute error percent V 
} MASS_CORRECT; ^ 

int NUMFAULTS; /* number of other detectable faults */ 

struct other_faultJnfb { 



* Each "other detectable fault has the following fields which 

* may be set as appropriate for the particular flowmeter application. 

* For example if a MAXimum parameter Is not relevant to detection 

* of a particular fault, it is set to +INFINITY. Similarly If a MINimum 

* parameter is not needed it is set to -INFINrTY. 



If the fault represented by the above signature is detected, then 

* the following fields define the actions to be taken with the pump 

* valve and panel display. 



float PUMPCHANGE; r adjustment to pump speed to correct */ 

float VALVECHANGE; /* adjustment to valve to correct V 

char *DISPLAY; r message to display */ 
} FAULT[MAXFAULTS]; 



V 

float MAXF; 
float MAXP; 
float MAXT; 
float MAXM; 
float MINF; 
float MINP; 
float MINT; 
float MINM; 
float MAXFSLOPE; 
float MAXPSLOPE; 
float MAXTSLOPE; 
float MAXMSLOPE; 
float MINFSLOPE; 
float MINPSLOPE; 
float MINTSLOPE; 
float MINMSLOPE; 
float MAXFCURV; 
float MAXPCURV; 
float MAXTCURV; 
float MAXMCURV; 
float MINFCURV; 
float MINPCURV; 
float MINTCURV; 
float MINMCURV; 

r 



r maximum mass flow slope */ 
r minimum frequency slope 7 
r minimum drive power slope V 
r minimum temperature slope */ 
r minimum mass flow slope 7 
r maximum frequency curvature 7 
r maximum drive power curvature 7 
r maximum temperature curvature 7 
r maximum mass flow curvature 7 
r minimum frequency curvature 7 
/* minimum drive power curvature 7 
/• minimum temperature curvature 7 
r minimum mass flow curvature 7 



/* maximum frequency 7 
/* maximum drive power 7 
r maximum temperature 7 
/* maximum measured mass flow 7 
r minimurn frequent 7 
r minimum drive power 7 
r minimum temperature 7 
r minimum measured mass flow 7 
/* maximum frequency slope 7 
/* maximum drive power slope 7 
r maximum temperature slope 7 



7 
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Fault Detection and Correction Methods - Overview: 

The methods of the present invention utilize the signature Information 
described above to detect and correct several fault conditions. FIG. 9 Is a flowchart 
of the methods of the present invention. The elements of FIG. 9 describe the 
5 methods of the present Invention at a high level. Later figures and associated 
descriptions below provide additional detail of the individual elements of FIG. 9. 

Element 900 operates to acquire samples of the current operating 
parameters of flowmeter 1 0 in preparation for analysis of the flowmeter operation. 
Element 900 is described below in further detail with respect to FIG. 10. 
10 Element 902 operates to store the samples acquired by element 900 in a 

circular buffer in RAM memory 201 . Element 902 is described below in further 
detail with respect to FIG. 1 1 . 

Element 904 operates to compute the slope and curvature of changes in the 
operating parameters sampled by operation of element 900 and stored in RAM 
15 memory 201 by operation of element 902. The most recent samples plus previous 
samples stored in the circular buffer in RAM memory 201 are used to determine the 
slope and curvature of changes in the operational parameters over time. Element 
904 is described below in further detail with respect to FIG. 12. 

Element 906 operates to determine whether a crack is forming in flow tubes 
20 12 and 14 of flowmeter 10. Element 906 analyzes the current sampled operating 
parameters and the slopes and curvatures computed in element 904 to detect the 
formation of a crack in flow tubes 12 and 14. If a crack is detected, element 906 
operates to stop the flow of material through flowmeter 10 and to relieve the fluid 
pressure within the flow tubes 12 and 14. Element 906 is described below in 
25 further detail with respect to FIG. 13. 

Element 908 operates to compute the corrected mass flow rate through 
flowmeter 10 by compensating the measured mass flow rate for the detected void 
fraction in the material flowing through flow tubes 12 and 14. If the conrected mass 
flow rate is found to be outside acceptable limits* element 908 operates to adjust 
30 the computed flow rate through flowmeter 10. Element 908 is described below in 
further detail with respect to FIG. 14. 

Element 910 operates to compare the sampled operating parameters and the 
computed slopes and curvatures of changes in the operating parameters against 
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corresponding signature information stored in ROM memory 202. if any fault 
signature information matches the con'esponding measured or computed 
information, then the connive actions described in ROM memory 202 along with 
the signature information are invoiced to attempt conrection of the problem and to 
5 notify an operator of the occurrence of the fault. Element 910 is described below 
in further detail with respect to FIG. 15. 

Element 912 operates to delay execution of the methods of the present 
invention until the next sampling period. This assures the periodic sampling and 
analysis of the operating parameters of flowmeter 1 0. Processing of the method 
1 0 continue at element 900. This method periodically repeats processing of elements 
900 through 912 to monitor operation of flowmeter 10 for purposes of detecting and 
correcting operating faults. 



Fault Detection and Correction Methods - Sampling: 

15 FIG. 10 is a flowchart depicting additional detail of the operation of element 

900 descrit)ed above. Elements of FIG. 10 operate to sample the operating 
parameters of flowmeter 10. Specifically, frequency, drive power, temperature, and 
measured mass flow are sampled by operation of the elements of FIG. 1 0 and 
stored in RAM memory 201 In variables named F, F. T, and fJl, respectively. 

20 Averaging of three consecutive readings helps to reduce the effects of noise in the 
analog signals received from mass flow instrumentation 24 over paths 160 through 
163. Variables stored in RAM memory 201 named SI, S2. and S3 are used to 
temporarily store three consecutive readings of each sampled operational 
parameter. 

25 Element 1000 operates to read three consecutive values of frequency from 

path 220 into variables SI, S2, and S3, respectively. 

Element 1002 computes the average frequency value by adding variables 

SI , S2. and S3 and dividing the sum by 3. The average frequency is then stored 

in the variable F for further processing by subsequent elements. 
30 Element 1 004 operates to read three consecutive values of drive power from 

path 221 into variables SI, S2, and S3, respectively. 
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Element 1006 computes the average drive power value by adding variables 
S1 , S2, and S3 and dividing the sum by 3. The average drive power is then stored 
in the variable F for further processing by subsequent elements. 

Element 1008 operates to read three consecutive values of temperature from 
5 path 222 into variables SI , S2, and S3, respectively. 

Element 1010 computes the average temperature value by adding variables 
SI , S2, and S3 and dividing the sum by 3. The average temperature is then stored 
in the variable T for further processing by subsequent elements. 

Element 1012 operates to read three consecutive values of measured mass 
10 flow from path 223 into variables S1, S2, and S3, respectively. 

Element 1014 computes the average measured mass flow value by adding 
variables SI, S2, and S3 and dividing the sum by 3. The average measured mass 
fbw is then stored in the variable fH for further processing by subsequent elements. 
This completes processing of the elements of FIG. 1 0 and therefor element 900 of 
15 FIG. 9. Processing continues with element 902 of FIG. 9. 



Fault Detection and Correction Methods - Storing: 

FIG. 1 1 is a flowchart depicting additional detail of the operation of element 
902 described above. Specifically, elements of FIG. 11 store the averaged 

20 operational parameter values computed by operation of element 900 of FIG. 9 in 
circular buffers in RAM memory 201 of FIG. 2. The values stored for the current 
sample period and previous sample stored in circular buffers are used in 
subsequent elements for further processing. 

Element 1 100 operates to store the averaged operating parameter values 

25 determined by operation of element 900 of FIG. 9: NEXT is a variable stored in 
RAM memory 201 and Is used as a circular index into parallel circular buffer arrays 
FARRAY, PARRAY, TARRAY, and MARRAY. The circular buffers are also stored 
in RAM memory 201. Specifically, F, P, T, and M are each stored in an element 
of FARRAY, PARRAY, TARRAY, and MARRAY. respectively, indexed by the NEXT 

30 variable current value. 

Element 1102 operates to increment the variable NEXT in a circular manner 
by applying a modulo function with MAXSAMPLES. MAXSAMPLES is the 
maximum number of samples to be stored in the circular buffer arrays discussed 
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above In element 1 100. This completes processing of the elements of FIG. 1 1 and 
therefor element 902 of FIG. 9. Processing continues with element 904 of FIG. 9. 

Fault Detection and Con-action Mathods - Slnp g and Curvature Compiitatinn,*;- 
5 FIG. 12 is a flowchart depicting additional detail of the operation of element 

904 described above. Elements of FIG. 12 operate to compute the slope, (first 
derivative with respect to time), and the curvature, (second derivative with respect 
to time), for changes in each stored average operational parameter. The variables 
used to store the computed slope and curvature of changes in each operation 
0 parameter are stored in RAM memory 201 . In addition, three temporary variables 
used to index into the circular buffers described above are stored in RAM memory 
201. 

Element 1200 computes indices into the circular buffer arrays FARRAY, 
PARRAY. TARRAY, and MARRAY for the past three average values stored In the 
5 circular buffere by operation of element 902 described above. Modulo arithmetic, 
as described above, is used to compute values for variables PREV1. PREV2, and 
PREV3. the indices into the circular buffers for the past three average values. 
PREV1, PREV2, and PREV3 are computed as the cun^nt value of NEXT minus 1, 
2, and 3 respectively. 

0 Element 1202 operates to compute the slope and curvature of changes in 

the stored averages of the frequency operational parameter. FSLOPE and FCURV 
are variables in RAM memory 201 of FIG. 2 in which the slope and curvature are 
stored. FSLOPE Is computed as the difference between the last two average 
values of frequency stored in FARRAYfPREVI] and FARRAY[PREV2]. The time 
base for the rate of change of the operational parameter is presumed to be one unit 
of time equal to the sampling period used to sample the operational parametere as 
discussed above in element 912 of FIG. 9. FCURV is the difference between the 
computed slopes of the past two time periods. Specifically. FCURV is equal to: 
(FARRAY[PREV1] - FARRAY[PREV2]) - (FARRAY[PREV2] - FARRAY[PREV3]) 
Element 1204 operates to compute the slope and curvature of changes in 
the stored averages of the drive power operational parameter. PSLOPE and 
PCURV are variables in RAM memory 201 of FIG. 2 in which the slope and 
curvature are stored. PSLOPE is computed as the difference between the last two 
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values of drive power stored in PARRAY[PREV11 arnJ PARRAY[PREV21. The time 
base for the rate of change of the operational parameter is presumed to be one unit 
of time equal to the sampling period used to sample the operational parameters as 
discussed above in element 912 of FIG. 9. PCURV is the difference between the 
5 computed slopes of the past two time periods. Specifically, PCURV is equal to: 
(PARRAYIPREV1] - PARRAYIPREV21) - (PARRAY[PREV2] - PARRAY[PREV3J) 
Element 1206 operates to compute the slope and curvature of changes in 
the stored averages of the temperature operational parameter. TSLOPE and 
TCURV are variables in RAM memory 201 of FIG. 2 in which the slope and 

10 curvature are stored. TSLOPE is computed as the difference between the last two 
values of temperature stored in TARRAY[PREV1) and TARRAY[PREV2I. The time 
t>ase for the rate of change of the operational parameter is presumed to be one unit 
of time equal to the sampling period used to sample the operational parameters as 
discussed above in element 912 of FIG. 9. TCURV is the difference between the 

1 5 computed slopes of the past two time periods. Specifically, TCURV Is equal to: 
(TARRAY[PREV1] - TARRAY[PREV2]) - (TARRAY[PREV21 - TARRAY[PREV3]) 
Element 1208 operates to compute the slope and curvature of changes In 
the stored averages of the measured mass flow rate operational parameter. 
MSLOPE and MCURV are variables in RAM memory 201 of FIG. 2 in which the 

20 slope and curvature are stored. MSLOPE is computed as the difference t)etween 
the last two values of measured mass flow rate stored in MARRAY[PREV1] and 
MARRAY[PREV2]. The time base for the rate of change of the operational 
parameter is presumed to be one unit of time equal to the sampling period used to 
sample the operational parameters as discussed above in element 912 of FIG. 9. 

25 MCURV is the difference between the computed slopes of the past two time 
periods. Specifically, MCURV is equal to: 

(MARRAY[PREV11 - MARRAY[PREV2]) - (MARRAYIPREV2] - MARRAYIPREV3]) 
This completes processing of the elements of FIG. 12 and therefor element 904 of 
FIG. 9. Processing continues with element 906 of FIG. 9. 

30 

Fault Detection and Correction Methods - Crack Detection: 

FIG. 13 is a flowchart depicting additional detail of the operation of element 
906 described above. Elements of FIG. 1 3 operate to determine whether a crack 
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is present in flow tubes 12 and 14 of flowmeter 10. A crack is detected by 
comparing the current sampled operational parameter values, computed slopes, 
and computed curvatures for frequency and drive power against signature 
information stored in ROM memory 202. REF.CRACK data structure described 
5 above contains the relevant limits for the particular flowmeter 10 in the particular 
application as shown in FIG. 1 A. Positively curved increasing drive power coupled 
with negatively curved decreasing frequency is the signature of a crack in a flow 
tube 12 or 14. 

Element 1300 operates to compare sampled drive power F, computed by 
0 operation of element 900, with REF.CRACK.MAXP discussed above. This 
determines whether drive power has increased beyond a threshold value such that 
further tests are required to determine that a crack exists in a flow tube 12 or 14. 
If the sampled driver power has not yet reached the corresponding signature 
threshold value, there is no crack in flow tubes 12 and 14 and processing of the 
5 elements of FIG. 12 is complete. If no crack is present, processing of element 906 
of FIG. 9 is complete and processing continues with element 906 of FIG. 9. 
Otherwise, further testing is require to determine whether a crad< is present in flow 
tubes 12 and 14 and processing continues with element 1302. 

Element 1302 operates to compare sampled frequency F, computed by 
operation of element 900, with REF.CRACK.MINF discussed above. This 
determines whether frequency has decreased beyond a threshold value such that 
further tests are required to determine that a crack exists in flow tube 12 and 14. 
If the sampled frequency has not yet reached the con-esponding signature threshold 
value there is no crack in flow tubes 12 and 14 and processing of the elements of 
FIG. 12 is complete. If no crack is present, processing of element 906 of FIG. 9 
is complete and processing continues with element 908 of FIG. 9. Otherwise, 
further testing is required to detemilne whether a crack is present in flow tubes 12 
and 14 and processing continues with element 1304. 

Element 1304 operates to compare the computed slope of changes in drive 
power PSLOPE, computed by operation of element 904. with 
REF.CRACKMAXPSLOPE discussed above. This determines whether drive power 
slope has increased beyond a threshold value such that further tests are required 
to determine that a crack exists in a flow tube 12 or 14. If the driver power slope 
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has not yet reached the corresponding signature threshold value there is no crack 
in flow tubes 12 and 14 and processing of the elements of FIG. 12 is complete, if 
no crack is present, processing of element 906 of FIG. 9 Is complete and 
processing continues with element 908 of FIG. 9. Othenfvise, further testing is 
5 require to determine whether a crack is present In flow tubes 12 and 14 and 
processing continues with element 1306. 

Element 1306 operates to compare the computed slope of changes in 
frequency FSLOPE, computed by operation of element 900, with 
REF.CRACK.MINFSLOPE discussed above. This detemiines whether frequency 

10 slope has decreased beyond a threshold value such that further tests are required 
to detenmine that a crack exists in a ftow tube 12 or 14. If the frequency slope has 
not yet reached the corresponding signature threshold value there is no crack in 
flow tubes 12 and 14 and processing of the elements of FIG. 12 is complete. If no 
crack Is present, processing of element 906 of FIG. i9 Is complete and processing 

15 continues with element 908 of FIG. 9. Othenvise, further testing is require to 
determine whether a crack is present in a flow tubes 12 or 14 and processing 
continues with element 1308. 

Element 1308 operates to compare the computed curvature of changes in 
drive power PCURV, computed by operation of element 904, with 

20 REF.CRACK.MAXPCURV discussed above. This determines whether drive power 
curvature has increased beyond a threshold value such that further tests are 
required to determine that a crack exists in flow tube 12 and 14. If the drive power 
curvature has not yet reached the con-esponding signature threshold value there is 
no crack in flow tubes 12 and 14 and processing of the elements of FIG. 12 Is 

25 complete. If no crack is present, processing of element 906 of FIG. 9 is complete 
and processing continues with element 908 of FIG. 9. Otherorise, further testing is 
require to determine whether a crack is present in flow tubes 12 and 14 and 
processing continues with element 1310. 

Element 1310 operates to compare the computed curvature of changes in 

30 frequency FCURV, computed by operation of element 900, with 
REF.CRACK.MINFCURV discussed above. This detemiines whether frequency 
curvature has decreased beyond a threshold value such that further tests are 
required to determine that a crack exists in flow tube 12 and 14. If the frequency 
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curvature has not yet reached the corresponding signature threshold value there is 
no crack In flow tubes 12 and 14 and processing of the elements of FIG. 12 is 
complete. If no crack is present, processing of element 906 of FIG. 9 is complete 
and processing continues with element 908 of FIG. 9. Otherwise, a crack is 
5 present in flow tubes 12 and 14 and processing continues with element 1312. 

Element 1312 operates to stop pump 180 by applying an appropriate signal 
to path 170. Stopping pump 180 will stop the flow of material through flow tubes 
12 and 14 to reduce the possibility of releasing material through the crack in flow 
tubes 12 and 14. 

! 0 Element 1314 operates to open valve 1 81 by applying an appropriate signal 

to path 171. Opening valve 181 reduces back pressure in flowmeter 10 to further 
reduce the possibility of releasing material through the crack in flow tubes 12 and 
14. 

Element 1316 operates to notify an operator of the existence of a crack in 
5 flow tubes 12 and 1 4 of flowmeter 1 0. An appropriate message is applied to path 
172 to be displayed on panel 182. This completes processing of the elements of 
FIG. 13 and therefor element 906 of FIG. 9. Processing continues with element 
908 of FIG. 9. 

Fault Detection and Con-action M ethods - Mass Flow Rata Correction- 

FIG. 14 is a flowchart depicting additional detail of the operation of element 
908 described above. Elements of FIG. 14 operate to compute the corrected mass 
flow rate through flowmeter 10. The corrected mass flow rate is adjusted by an 
en^r factor (error percentage) determined from the measured mass flow rate and 
a computed void fraction of the material flowing through flow tubes 12 and 14. Void 
fraction is determined as a function of sampled frequency and sampled drive power 
as depicted in FIG. 6 and as discussed above. Error percentage is determined as 
a function of the computed void fraction and the measured mass flow provided by 
mass flow instmmentation 24 as discussed above with respect to FIG. 8. The error 
percentage is then applied to the measured mass flow rate to determine the 
corrected mass flow rate. The Information shown in FIGS. 6 and 8 is encoded in 
ROM memory 202 as discussed above for use in determining void fraction and 
error percentage to correct the measured mass flow rate. 
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REF.MASS_CORRECT.VOID_FRACTION is the encoding of the Information stored 
in ROM memory 202 and shown in FIG. 6 used for computing void fraction as a 
function of measured frequency and drive power. 
REF.MASS_CORRECT.ERROR_PERCENT is the encoding of the information 
5 stored in ROM memory 202 and shown In FIG. 8 used for computing error percent 
as a function of measured mass flow rate and computed void fraction. 

Use of corrected mass flow rate enables more precise measurement under 
a wider range of process conditions. 

Element 1400 operates to detennine the void fraction for the material flowing 

10 through flowmeter 10. Void fraction is determined as a function of F, F, FSLOPE, 
PSLOPE, FCURV. and PCURV. all detemiined by operation of elements 900 and 
904 discussed above, used in conjunction with the Information depicted In FIG. 6 
and encoded in ROM memory 202 as REF.MASS_CORRECT.VOID_FRACTION. 
As discussed above, the information In ROM memory 202 may be encoded In 

15 tabular form or as coefficients of a polynomial equation used to determine void 
fraction. Tabular fbmi is preferred for low cost microprocessors such as CPU 200 
of FIG. 2 with limited computational capabilities. The current values of F, P, 
FSLOPE, and PSLOPE are used to interpolate in the table of entries stored at 
REF.MASS_CORRECT.VOID_FRACTION in ROM memory 202 to determine the 

20 current void fraction of the material flowing within flowmeter 10. 

Element 1402 operates to determine the error percent to be applied to 
correct the measured mass flow rate of flowmeter 10. En-or percent is determined 
as a function of M, determined by operation of element 900 discussed above, and 
the void fraction, computed above by element 1400, used in conjunction with the 

25 infonmation depicted in FIG. 8 and encoded in ROM memory 202 as 
REF.MASS_CORRECT-ERROR_PERCENT. As discussed above, the infomiation 
in ROM memory 202 may be encoded In tabular form or as coefficients of a 
polynomial equation used to determine void fraction. Tabular forni is preferred for 
low cost microprocessors such as CPU 200 of FIG. 2 with limited computational 

30 capabilities. The current values of Kland the void fraction determined by operation 
of element 1400 above are used to interpolate in the table of entries stored at 
REF.MASS_CORRECT.ERROR„PERCENT in ROM memory 202 to determine the 
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current error percentage for correcting the measured mass flow rate of material 
flowing within flowmeterlO. 

Element 1404 operates to determine the corrected mass flow rate from ¥A 
and the computed error percentage determined by operation of element 1402 
discussed above. The emor percentage is applied to the measured mass flow rate 
to correct for the effects of void fraction on the measured mass flow rate. 

Element 1406 operates to notify the operator of the corrected mass flow rate 
by displaying a message on panel 182 sent over path 172. 

Element 1408 compares the corrected mass flow rate, determined by 
operation of element 1406, with a nominal threshold minimum flow rate stored in 
ROM memory 202 at REF.MINM, as discussed above. If the corrected mass flow 
rate is below the threshold value, processing continues with element 1410. 
Othenwise, processing continues with element 1412. 

Element 1410 operates to increase the mass flow rate by increasing the 
speed of pump 180 and by opening valve 181. CPU 200 of FIG. 2 controls pump 
180 and valve 181 by applying signals to paths 170 and 171 respectively. This 
completes processing of the elements of FIG. 14 and therefor element 908 of FIG. 
9. Processing continues with element 910 of FIG. 9. 

Element 1412 compares the corrected mass flow rate, determined by 
operation of element 1406. with a nominal threshold maximum flow rate stored in 
ROM memory 202 at REF.MAXM, as discussed above. If the corrected mass flow 
rate is above the threshold value, processing continues with element 1414. 
Othenwise, this completes processing of the elements of FIG. 14, and therefor 
element 908 FIG. 9, and processing continues with element 910 of FIG. 9. 

Element 1414 operates to decrease the mass flow rate by decreasing the 
speed of pump 180 and closing valve 181. CPU 200 of FIG. 2 controls pump 180 
and valve 181 by applying signals to paths 170 and 171 respectively. This 
completes processing of the elements of FIG. 14 and therefor element 908 of FIG. 
9. Processing continues with element 910 of FIG. 9. 



Fault Detection and Correction Met hods - Other Fault Detection and Cormntinn- 

FIG. 15 is a flowchart depicting additional detail of the operation of element 
910 described above. Elements of FIG. 15 operate to detect other operating fault 



wo 96/05484 PCT/US95/09324 

-31- 

condltions in flowmeter 10. As discussed al30ve, signature information for each 
detectable fault condition is stored in ROM memory 202 encoded as high and low 
threshold values for each sampled operating parameter as well as high and low 
threshold values for the associated computed slope and curvature of changes in 
5 each corresponding operational parameter. The number of other fault conditions 
detectable by the methods of the present invention are stored in ROM memory 202 
in REF.NUMFAULTS as discussed above. The signature information for all 
detectable fault conditions Is stored in the array of staictures REF.FAULT[ ] as 
discussed above. 

1 0 One skilled in the art will recognize that detection of a crack in flow tubes 1 2 

and 14 of flowmeter 10 as discussed above with reference to FIG. 13, is an 
example of a fault detectable under operation of the general fault detection methods 
described herein with respect to FIG. 15. Crack detection is separately addressed 
above with reference to FIG. 13 to help clarify the precise method of detecting a 

1 5 crack under the methods of the present invention. The threshold values encoded 
in REF.CRACK may be similarly encoded as one of the faults detectable by 
operation of the methods of FIG. 15 and encoded in ROM memory 202 within 
REF.FAULT. 

Element 1500 operates to initialize to zero a loop counting variable, named 
20 I, stored in RAM memory 201. I is used in subsequent elements to detemiine when 
processing of all detectable fault conditions is completed. 

Element 1502 operates to compare each sampled operating parameter to its 
corresponding high and low threshold value encoded and stored in ROM memory 
202. Specifically, M Is compared to REF.FAULT[I].MINM and to 
25 REF.FAULT[I].MAXM. In like manner, P, F, and T are each compared to their 
respective minimum and maximum threshold values, namely: REF.FAULT[I].MINP, 
REF.FAULT[I].MAXP. REF.FAULT[I].MINF. REF.FAULT[I].MAXF, 
REF.FAULTII].MINT, and REF.FAULT[I].MAXT. If any comparison indicates that 
any one of the sampled operating parameters is outside the specified range, the 
30 corresponding fault has occurred and processing continues with element 1512. 
Otherwise, processing continues with element 1504. 

Element 1504 operates to compare the computed slope of the change in 
each sampled operating parameter to its corresponding high and low threshold 
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value encoded and stored in ROM memory 202. Specifically, MSLOPE is 
compared to REF.FAULTIIJ.MINMSLOPE and to REF.FAULT[l].MAXMSLOPE. In 
like manner, PSLOPE, FSLOPE, and TSLOPE are each compared to their 
respective minimum and maximum threshold values, namely: 
5 REF.FAULT[l].MINPSLOPE. REF.FAULT[l].MAXPSLOPE, 
RE F.FAULT[I].M IN FSLOPE, REF.FAULT[l].MAXFSLOPE. 
REF.FAULTIIJ.MINTSLOPE, and REF.FAULT[I].MAXTSLOPE. If any comparison 
indicates that any one of computed slopes is outside the specified range, the 
corresponding fault has occurred and processing continues with element 1512. 

10 OthenAfise, processing continues with element 1506. 

Element 1506 operates to compare the computed curvature of the change 
in each sampled operating parameter to Its corresponding high and low threshold 
value encoded and stored in ROM memory 202. Specifically, MCURV is compared 
to REF.FAULT[I1.MINMCURV and to REF.FAULT[I].MAXMCURV. In like manner, 

1 5 PCURV, FCURV, and TCURV are each compared to their respective minimum and 
maximum threshoW values, namely: REF.FAULT[I1.MINPCURV, 
REF.FAULT[I].MAXPCURV. REF.FAULT[I].MINFCURV, REF.FAULT[I].MAXFCURV, 
REF.FAULTII].MINTCURV, and REF.FAULT[I].MAXTCURV. If any comparison 
indicates that any one of the computed curvatures is outside the specified range, 

20 the corresponding fault has oocun-ed and processing continues with element 1512. 
Otherwise, processing continues with element 1508. 

Element 1508 operates to increment loop counter I by one to prepare for 
testing the next fault condition descril>ed in the REF.FAULTI ] anay of encoded 
signature information. 

25 Element 1510 operates to compare loop counter I to REF.NUMFAULTS. If 

I has not yet reached the number of faults encoded in the REF.FAULTI ] array, then 
processing continues with element 1502 to process another fault signature encoded 
in ROM memory 202. Othenwse, this completes processing of the elements of FIG. 
15 and therefor element 910 of FIG. 9. Processing continues with element 912 of 

30 FIG. 9. 

Element 1512 operates to change the speed of pump 180 according to the 
percentage value in REF.FAULT[I].PUMPCHANGE. The speed of pump 180 is 
changed by CPU 200 applying an appropriate signal on path 170. The change in 



wo 96/05484 PCT/US95/09324 

-33- 

speed of pump 180 is intended to alter the mass flow rate of material In flowmeter 
10 in order to correct the detected fault condition encoded In REF.FAULT[I]. 

Element 1514 operates to change the position of valve 181 according to the 
percentage value in REF.FAULTIi].VALVECHANGE. The position of valve 181 Is 
5 changed by CPU 200 applying an appropriate signal on path 171. The change in 
the position of valve 181 is intended to alter the mass flow rate of material in 
flowmeter 10 in order to correct the detected fault condition encoded in 
REF.FAULT[I]. 

Element 1516 operates to display a message on panel 181 indicating the 

10 occurrence of the detected fault condition. The message to be displayed is 
retrieved from the value in REF.FAULTIIJ.DISPLAY. The message is displayed on 
panel 182 by CPU 200 applying an appropriate signals on path 172. This 
completes processing of the elements of FIG. 15 and therefor element 910 of FIG. 
9. Processing continues with element 912 of FIG. .9. 

15 It Is expressly understood that the claimed invention is not to be limited to 

the description of the preferred embodiment but encompasses other nKxJifications 
and alterations within the scope and spirit of the Inventive concept. In particular, 
the methods and apparatus of the present invention may be applied to Coriolis 
effect mass flowmeters with a variety of tube shapes including so called "U" shaped 

20 tubes, straight Xubes. and others. In addition, the information and data shown on 
the graphs of FIGS. 3 through 8 represent measurements made on a plurality of 
flowmeters, at different times, and under different operating conditions. The data 
from one graph of FIGS. 3 through 8 therefore cannot be meaningfully compared 
to the data of another graph. There are many variables involved in detennining 

25 flowmeter data output. Among these variables are: thickness of the flow tube walls, 
material temperature, viscosity, and flowrate: elastic properties of the material 
comprising the flow tubes; etc. Therefore, the data shown on each of these figures 
is only illustrative and should be used only to portray the relationships between the 
various parameters shown on the same graph. Because of the many variables 

30 involved, it cannot t>e assumed that the numerical values shown will be readily 
reproduced by others. 
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WE CLAIM: 

1. A method for detecting changes in the operating condition of an 
operating Coriolis effect mass flowmeter (10), CHARACTERIZED IN THAT said 
method comprises the steps of: 

detennining at least one value indicative of present operation of said 
flowmeter selected from the group consisting of present operating parameters, 
slopes of the changes in said present operating parameters, and curvatures of said 
slopes (900, 904); 

determining that said at least one value has deviated from a nominal 
value by a predetermined amount (906, 908, 910); 

generating an output signal indicating the presence of said change in 
said at least one value in response to said determination that said at least one 
value has deviated from a nominal value by said predetermined amount (906, 908. 
910); and 

applying said output signal to utilization means for purposes including 
the correction of said abnormal condition in said flowmeter (906, 908, 910). 

2. The method of daim 1 wherein said step of applying said output signal 
to said utilization means comprises controlling a pump (180) associated with said 
flowmeter to alter the flow rate of material flowing through said flowmeter. 

3. The method of claim 1 wherein said step of applying said output signal 
to said utilization means comprises controlling a valve (181) assodated with said 
flowmeter to alter the flow rate of material flowing through said flowmeter. 

4. The method of claim 1 wherein: 

said nominal value of said stope of said vibration frequency is greater 
than or equal to ad^out zero (1302); 

said nominal value of said slope of said drive power is less than or 
equal to about zero (1300); 

said nominal value of said cun/ature of said vibration frequency is 
greater than or equal to about zero (1310); 
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said nominal value of said curvature of said drive power is less than 
or equal to about zero (1308); and 

said output signal indicates the presence of a crack in said flow tubes 
of said flowmeter (1312, 1314, 1316). 

5- The method of claim 4 wherein said step of applying said output signal 
to said utilization means comprises controlling a pump (180) associated with said 
flowmeter to stop material from flowing through said flowmeter (1312). 

6. The method of claim 4 where said step of applying said output signal 
to said utilization means comprises controlling a valve (181) associated with said 
flowmeter to stop material from flowing through said flowmeter (1314). 

7. The method of claim 1 wherein: 

said nominal value of said slope of said vibration frequency is greater 
than or equal to zero (1504); 

said nominal value of said slope of said drive power is zero (1504); 

said nominal value of said curvature of said vibration frequency is less 
than or equal to zero (1506); and 

said output signal indicates the presence of an unacceptable density 
change In material flowing in said flow tubes of said flowmeter (1512, 1514, 1516). 

8. The method of claim 1 wherein: 

said nominal value of said slope of said vibration frequency is less 
than or equal to zero (1504); 

said nominal value of said slope of said drive power is less than or 
equal to zero (1504); 

said nominal value of said curvature of said vibration frequency is zero 

(1506); 

said nominal value of said curvature of said drive power is less than 
or equal to zero (1506); and 

said output signal indicates the presence of an unacceptable void 
fraction in material flowing in said flow tubes of said flowmeter (1512. 1514, 1516). 
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9. The method of claim 1 wherein: 

said nominal value of said slope of said vibration frequency is greater 
than or equal to zero (1504); 

said nominal value of said slope of said drive power is less than or 
equal to zero (1504); 

said nominal value of said curvature of said vibration frequency is zero 

(1506); 

said nominal value of said curvature of said drive power is zero 

(1506); and 

said output signal indicates the presence of an unacceptable mass 
fraction in material flowing in said flow tubes of said flowmeter (1512, 1514, 1516). 

10. The method of claim 1 further comprising the steps of: 
determining a void fraction value of material flowing through said flow 

tut>es as a function of said at least one operating parameter (1400); 

determining an error correction factor as a function of said present 
measured mass flow rate and said void fraction value (1402); 

determining a corrected mass flow rate using said present measured 
mass flow rate and said error correction factor (1404); and 

applying an output signal proportional to the difference between said 
con-ected mass flow rate and said present measured mass flow rate to utilization 
means for altering the rate of flow of material flowing through said flowmeter (1406, 
1408, 1410, 1412, 1414). 

1 1 . The method of claim 10 wherein said step of detennlning said void 
fraction value further comprises interpolating within a table of values to determine 
said void fraction value as a function of said present vibration frequency. 

12. The method of daim 10 wherein said step of detemiining said void 
fraction value further comprises interpolating within a table of values to determine 
said void fraction value as a function of said present drive power. 
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is. The method of claim 10 wherein said step of determining said error 
correction factor comprises interpolating within a table to determine said error 
correction factor as a function of said present measured mass flow rate and of said 
void fraction value. 

14. An apparatus for detecting a fault condition in a Coriolis effect mass 
flowmeter (10), CHARACTERIZED IN THAT said apparatus further comprises: 

input means connected to said flowmeter for receiving information 
regarding present operating parameters of said flowmeter (160-163. 203-206); 
5 memory means for storing fault signature information regarding at 

least one fault condition for said present operating parameters (201 , 202); 

output means for indicating the occurrence of said fault condition in 
operation of said flowmeter (180, 181, 182); and 

detection means connected to said input means, connected to said 
10 output means, and connected to said memory means, for comparing said present 
operating parameters to said fault signature information to detect the occurrence for 
said fault condition and to apply a signal to said output means indicating the 
occurrence of said fault condition (200). 
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